Determination of the total width of the rj meson 
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Taking advantage of both the low-emittance proton beam of the cooler synchrotron COSY and 
the high momentum precision of the COSY-11 detector system, the mass distribution of the 77' 
meson was measured with a resolution of 0.33 MeV/c 2 (FWHM), improving the experimental mass 
resolution by almost an order of magnitude with respect to previous results. Based on the sample 
of more than 2300 reconstructed pp — > ppn' events the total width of the rj' meson was determined 
to be F v , = 0.226 ± 0.017(stat.) ± 0.014(syst.) MeV/c 2 . 

PACS numbers: 13.60.Le, 14.40.Bc, 14.70.Dj 
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In this Letter, we report on the measurement of the 
mass distribution of the 77' meson carried out with a res- 
olution of a fraction of MeV/c 2 . This accuracy was ob- 
tained by using the low-emittance proton-beam of the 
cooler synchrotron COSY [l[ and the high momentum 
resolution of the COSY-11 detector system [1, and it 
is nearly an order of magnitude more precise than previ- 
ous results. 

In the latest review by the Particle Data Group 
(PDG) i, two values for the total width of the 77' meson 
are given. One of these values, (0.30±0.09) MeV/c 2 , re- 



sults from the average of two measurements [5|,[6|], though 
only in one of these experiments was T„/ extracted di- 
rectly based on the mass distribution [5J. The second 
value (0.205±0.015) MeV/c 2 , recommended by the PDG, 
is determined by fit to altogether 51 measurements of par- 
tial widths, branching ratios, and combinations of parti- 
cle widths obtained from integrated cross sections Q. 
The result of the fit is strongly correlated with the value 
of the partial width r(7y' — > 77), which causes serious dif- 
ficulties when the total and the partial width have to be 
used at the same time, like e.g. in studies of the gluonium 
content of the 77' meson 

The partial width of the 77' — > 77 channel can be ex- 
tracted from the e + e~ — >• e + e~?7' cross sections without 
knowledge of the 0- 11|, yet its derivation is model- 
dependent due to the need to incorporate a form fac- 
tor which describes the spatial distribution of the elec- 
tric charge in the 77' meson. For the derivation of par- 
tial widths of all other decay channels, the knowledge of 
r, ( ' is mandatory. At present it is the inaccuracy of T v i 
which limits investigations of many interesting physics 
issues, such as, for example, the quark mass difference 
mrid — rn u 12? , 13 1, isospin breaking in QCD [12, 14 1. 



or the box anomaly of QCD [15j. This is because the 
branching ratios of the 77' meson decay channels are typi- 
cally known with a relative precision of more than an or- 
der of magnitude better than the present accuracy with 



which r i; / is extracted 



The signal of the 7/ meson production observed in 
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previous experiments (ol. lil. Il6421l| with mass resolutions 
poorer than ~1 MeV/c 2 do not a priori exclude the pos- 
sibility that some structure in the mass distribution of 
the 77' meson would be visible at higher precision. Ex- 
tractions of the 77' width (r, ( / ) were performed under the 
assumption that the 77' meson is a single state. This, how- 
ever, must not necessarily be the case [22, 23 1 if there is a 
significant glue contribution in the wave function of this 
meson j24[. The precision achieved with the COSY-11 fa- 
cility enabled us for the first time to determine the mass 
spectrum of the 77' meson with a resolution comparable 
to its total width of -0.2 MeV/c 2 

The experiment, reported in this Letter, was per- 
formed in the Research Centre Julich. The value of 
r^/ was established directly from the measurement of 
the mass distribution of the 77' meson, produced via the 
pp — > ppij reaction. The momentum of the COSY beam 
and the dedicated zero degree COSY-11 facility enabled 
the measurement at an excess energy of only a fraction 
of an MeV above the kinematic threshold for the 77' me- 
son production. This was the most decisive factor in 
minimizing uncertainties of the missing-mass determina- 
tion, since at threshold the partial derivative of the miss- 
ing mass with respect to the outgoing proton momentum 
tends to zero. In addition, close to threshold the signal- 
to-background ratio increases due to the more rapid re- 
duction of the phase space for multimeson production 
than for the 77'. 
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In order to control systematic uncertainties, the mea- 
surement was carried out at five different beam mo- 
menta, which were nominally 3211, 3213, 3214, 3218, 
and 3224 MeV/c. The cooled beam of protons [H? ] 
circulated in the ring of the cooler synchrotron COSY 
through a stream of the hydrogen cluster target [2(| . In 
the magnetic field of the COSY dipole, the final state pro- 
tons from the pp — > pprj' reaction were bent more than 
the beam protons and were measured by means of the 
COSY-11 detector shown schematically in Fig. [TJ The 
momentum vectors of the outgoing protons were recon- 
structed based on the bending of their trajectories in the 
magnetic field between the center of the reaction region 
and the tracks measured in the drift chambers (Dl and 
D2). In addition, the velocities of the two protons were 
determined from their time of flight measured between 
the scintillator detectors SI and S3. The independent 
determination of momentum and velocity enables parti- 
cle identification to be made via its invariant mass. Since 
the momentum is reconstructed more precisely than the 
velocity, after the identification, the energy of the parti- 
cle is derived from its known mass and momentum. The 
standard technique for monitoring the beam momentum 
at the COSY accelerator is through the measurement of 
the frequency distribution of the circulating beam. Such 
a distribution can be transformed to the momentum coor- 
dinate by using the values of the accelerator settings [25| . 
As an example, a spectrum for the lowest beam energy 
used in the experiment is presented in the left corner 
of Fig. Q] The beam momentum distribution is smooth 
and its spread is equal to 2.5 MeV/c (FWHM). How- 
ever, due to the position of the COSY-11 target system 
in a bending section of the COSY ring in a dispersive 
region, the effective spread of the beam (the momentum 
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Figure 1: Schematic view of the COSY-11 detector setup 
(top view). SI, S2, S3 and S4 denote scintillator detectors, Dl 
and D2 indicate drift chambers and Si stands for the silicon- 
pad detector. Left corner: Momentum spectrum for the mea- 
surement with the nominal beam momentum of 3211 MeV/c. 
As an example, the effective spread of the beam momentum 
due to the dispersion is shown as for a target width of 1 mm 
(dashed line) and 1 cm (solid line). 



range seen by the target) is smaller [28[. The 7/ meson 
was not registered but instead it was identified by using 
the missing-mass technique. The precision of the deter- 
mination of the size and position of the target stream 
influences the accuracy of the reconstruction of the mo- 
mentum of the outgoing particles and the accuracy of the 
determination of the momentum spread of beam protons 
interacting with the target. As a compromise between 
accuracy and statistics, the transverse size of the target 
stream was reduced to 0.9 mm which is significantly less 
than the horizontal spread of the COSY beam. There- 
fore, the momentum spread of the interacting protons is 
defined by the momentum dispersion at the target region 
and by the size of the target stream. The size and posi- 
tion of the target stream, being crucial for the analysis, 
were monitored by two independent methods. The first 
was based on the measurement of the momentum dis- 
tribution of elastically scattered protons |27j |. while the 
second was a direct measurement of the target geome- 
try by mechanically scanning the target stream position 
above and below the target area from time to time. A 
diagnostic unit with several wires was rotated through 
the target stream and the pressure in the cluster beam 
dump was measured as a function of the wire position. 
When parts of the target stream are blocked by a wire, 
the pressure decreases proportional to the blocked area. 
Therefore variations of the pressure allowed the monitor- 
ing of the size and alignment of the target stream dur- 
ing the experiment. The results of the two methods are 
in good agreement and the achieved precision is ±0.05 
and ±0.01 mm for the size and alignment, respectively. 
The momentum distributions of the elastically scattered 
protons were used not only for monitoring the relative 
geometrical settings of the target, but also for the exact 
positions of the dipole field and the drift chambers [28| . 

As the next step of the analysis, the missing- mass spec- 
tra were determined in order (i) to distinguish between 
the signal and background, (ii) to evaluate the absolute 
beam momenta, and finally (hi) to extract the width of 
the 77' meson. The spectra for the highest and lowest 
excess energies are shown in Fig [2j It is important to 
stress that the background distribution is smooth in the 
whole range studied, and that the signal from rf me- 
son production shows up clearly. The figure illustrates 
that the spectrum at one energy can be used as a good 
estimate of the background to the spectra at the other 
energies. The method for the background subtraction is 
based on the observation that the shape of the multip- 
ion mass distribution does not change when the excess 
energy for the pp — >• pprj' reaction varies by a few MeV, 
which is small compared to the total available energy of 
about 500 MeV [l6| . The systematic error in the changes 
of the shape due to the method applied was estimated to 
be less than 1% even for shifts many times larger than 
the energy range relevant in these measurements [29| . In 
order to decrease the influence of the statistical fluctua- 
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tion, the background for a given energy was taken from 
a second-order polynomial fit to the data at a different 
energy, which was shifted and normalized to the data of 
interest [II]. 

Owing to the large statistics of the momentum distri- 
butions for elastically scattered protons, relative differ- 
ences between excess energies were determined with a 
negligible statistical error from the sizes of the kinematic 
ellipses [281 ], Next, the absolute values of the excess ener- 
gies were derived by comparing the position of the mean 
of the missing- mass peak for data closest to the threshold 
with the empirical value of the mass of the rj' meson Q ■ 
The "true" values of the beam momenta thus determined 
are 3210.7, 3212.6, 3213.5, 3217.2, and 3223.4 MeV/c, 
corresponding to excess energies of 0.8, 1.4, 1.7, 2.8, 
and 4.8 MeV, respectively. The accuracy of the beam 
momentum determination amounts to ±0.2 MeV/c and 
is predominantly due to the uncertainty of the rj' mass 
(957.78 ± 0.06) MeV/c 2 0. The systematically lower 
values of the true beam momenta of about 0.5 MeV/c 
are consistent with previous experience at COSY where 
the real beam momentum was always smaller than the 
nominal value (30I ? ] . 

In order to derive the value of the rj width, the exper- 
imental missing-mass spectra were compared with dis- 
tributions simulated with different values of IV . In the 
simulations based on the GEANT3 packages [31| , the re- 
sponse of the COSY- 11 detector system to the pp — > pprj' 
reaction was generated, taking into account the geome- 
try and material composition, as well as relevant reso- 
lutions of the COSY-11 detector components, including 
the size of the target stream, the spatial and momentum 
spread of the beam, and also all known physical processes 
such as multiple scattering and nuclear reactions. In the 
simulations of the mass distribution, a Breit-Wigner for- 
mula for the 77' meson was used. Afterwards the gener- 
ated events were analyzed in the same way as the ex- 
perimental data, and sets of missing-mass spectra were 
reconstructed for the values of T n i ranging from 0.14 to 
0.38 MeV/c 2 . Finally, the sum of the experimental back- 
ground and the Monte Carlo missing-mass spectra for 
the pp — > ppr/' reaction was fitted to the experimental 
data. The normalization factor of the r/ signal was the 
only free parameter in this fit. The result of the fit is 
shown by solid lines in Fig. [3] The decrease of the width 
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Figure 2: Missing-mass spectra for the pp — > ppX reac- 
tion determined at beam momenta of 3211 (open points) and 
3224 MeV/c (filled points). The filled points were shifted to 
the kinematic limit and normalized to the open points. 



of the missing-mass with decreasing excess energy is a 
kincmatical effect reflecting the propagation of errors of 
momenta involved in the missing mass calculations [321 ]. 
The simulations reproduce very well the change of the 
signal width with excess energy and thus validate the 
correctness of the established detector and target char- 
acteristics. The lower-right panel of Fig. [3] presents the 
dependence of the x 2 011 the T v i value. The minimum 
of x 2 is at IV = 0.226 MeV/c 2 and the la statistical 
error is equal to ±0.017 MeV/c 2 . 

The systematic error was estimated by studying the 
sensitivity of the result to the variation of parameters de- 
scribing the experimental conditions in the analysis and 
in the simulation [28j . The contributions to the system- 
atic error are (i) the target position (±0.006 MeV/c 2 ) and 
size (±0.002 MeV/c 2 ), (ii) the position and orientation 
of the drift chambers (±0.001 MeV/c 2 ), (iii) the map of 
the magnetic field (±0.007 MeV/c 2 ), and (iv) the abso- 
lute beam momentum determination (±0.003 MeV/c 2 ). 
These values were estimated as the difference between the 
derived result of the T v > and the T 7 y values established 
by changing in the analysis and simulations a particular 
parameter by its error. The systematic error due to the 
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Figure 3: The missing-mass spectra for the pp — > ppX reac- 
tion. The r\ meson signal is clearly visible. The experimental 
data are presented as points, while in each plot the line corre- 
sponds to the sum of the Monte Carlo generated signal for the 
pp — t ppvl reaction with T v i = 0.226 MeV/c 2 and the back- 
ground obtained from another energy. The plot at the bottom 
right of the figure presents x 2 as a function of the T v i. The 
minimum value of the x 2 divided by the number of degrees of 
freedom amounts to 0.96. 
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Figure 4: Comparison of available values of T v i . 

method of the background subtraction (±0.006 MeV/c 2 ) 
was established as the maximum difference between T v i 
values determined when using experimental background 
shapes from different energies. The uncertainty due 
to the bin width (±0.004 MeV/c 2 ) was estimated by 
changing the width of bins in the range from 0.1 to 
0.04 MeV/c 2 . Furthermore, the sensitivity of the result 
to the range of the missing-mass values used for the fit 
(±0.005 MeV/c 2 ) was estimated by enlarging the mass 
range by seven bins on each side of the peak. The inaccu- 
racy due to the model applied in the simulations for the 
proton-proton final state interaction (±0.003 McV/c 2 ) 
was estimated conservatively as a differences in results 
determined when using parameterization of the proton- 
proton S-wave interaction [33l |3~H and when neglecting 
the final-state interactions. Finally, the total system- 
atic error was estimated as the quadratic sum of the 
nine independent contributions mentioned above and is 
0.014 MeV/c 2 . Our final result is compared with earlier 
width determinations in Fig. 2] 

In summary, the mass distribution of the ?/ meson 
has been measured with an experimental resolution of 
FWHM = 0.33 McV/c 2 . The rf meson was created in 
the pp — > pprj' reaction close to the kinematic threshold 
by using the low-emittance proton beam of the cooler 
synchrotron COSY incident on a stream of hydrogen clus- 
ters. The outgoing protons were detected by using the 
COSY- 11 facility. The total width of the rj meson was 
extracted from the missing-mass spectra and amounts to 
T v > = 0.226 ± 0.017(stat.) ± 0.014(syst.) MeV/c 2 . The 
result does not depend on knowing any of the branch- 
ing ratios or partial decay widths. The extracted T v > 
value is in agreement with both previous direct deter- 
minations of this value (r„/ = 0.28 ± 0.10 McV/c 2 @ 
and = 0.40 ± 0.22 MeV/c 2 0). The achieved accu- 
racy is similar to that obtained by the PDG from a fit to 
51 measurements of branching ratios and cross sections 
(IV = 0.204 ± 0.015 MeV/c 2 ) 0. 

Finally, it is worth noting that the achieved mass res- 
olution is of the same order as the total width of the 77' 
meson itself, thereby excluding the possibility of a sub- 
structure in the rf signal at this level. 
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